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Abstract. In the framework of the A;7^-factorization approach, we analyse the inclusive and 
inelastic production of J/V' particles in deep inelastic ep scattering. We take into account 
both colour-singlet and colour-octet production channels. We inspect the sensitivity of theo- 
retical predictions to the choice of model parameters. Our theoretical results agree reasonably 
well with recent experimental data collected by the collaboration HI at HERA. 

1 Introduction 

Over the last decade, investigation of the J/ip production mechanisms in hadron- 
hadron and lepton-hadron collisions continues to attract significant attention from both 
theoretical and experimental sides. The puzzling history of J/ip traces back to the 
early 1990s, when the measurements [l]-[4] of the J/tp and T hadropro duct ion cross 
sections revealed a more than one order-of-magnitude discrepancy with theoretical ex- 
pectations [5]- [7]. This fact has induced extensive theoretical activity. In particular, 
it led to the introduction of a new production mechanism, the so called colour-octet 
model [8]- [16]. Since then, the colour-octet model has been believed to give the most 
likely explanation of the quarkonium production phenomena, although there are also 
some indications that it is not working well [17]. One of the problems is connected 
with the photoproduction data [18], [19] where the contribution from the colour-octet 
mechanism is unnecessary or even unwanted [20]- [22] as the experimental results can be 
described within the colour- singlet model alone (if the next-to-leading-order contribu- 
tions are taken into account) [23]. Another difficulty refers to the J/tp spin alignment. 
If, as expected, the dominant contribution comes from the gluon fragmentation into an 
octet cc pair, the mesons must have strong transverse polarization [22]- [26]. This is 
in disagreement with the data [27], [28], [29], [30] which point to unpolarized or even 
longitudinally polarized mesons. 

A different strategy is represented by the /cT^-factorization approach [31]- [34] (see 
also recent review [35] and references therein). In this approach, one focuses on the re- 
summation of "small x" logarithms (i.e., the terms [ln(/i^/A^) Og]^, [ln(/i^/A^) ln(l/a;) a^] 
and [ln(l/a;) ««]"') to all orders in n and takes into account the effects of finite transverse 
momenta of partons. The resummation results in the "unintegrated" parton distribu- 
tions which generalize the QCD factorization beyond the coUinear approximation. 

The effects of the initial gluon transverse momentum and gluon off-shellness were 
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shown to have an important impact on the J /ip production properties in the photon- 
gluon fusion colour-singlet model [36], [37]. Several attempts have been made in the 
literature to incorporate the colour-octet model within the /cy-factorization scheme 
[38], [39]. An extensive analysis of the production of J/V', Xc and T mesons in pp 
collisions at the Fermilab Tevatron has been recently presented in Ref. [40]. Our 
paper is a logical continuation of this line. Here, we concentrate on the inclusive and 
inelastic production of J/ if] particles by virtual photons in deep inelastic ep scattering 
at HERA. The calculations are based on the /cT-factorization approach and on the 
nonrelativistic QCD formalism where we take into account both colour-singlet and 
colour-octet channels. 

The outline of the article is as follows. First, in Sec. 2, we briefly recall the basic 
principles of the colour-singlet and colour-octet models and explain their extension to 
the /cT-factorization approach. The necessary technical details are described in Sec. 3. 
The numerical results followed by a discussion are displayed in Sec. 4. Finally, the 
concluding remarks are collected in Sec. 5. 

2 Theoretical framework 

In the framework of the colour-singlet approach [5]-[7], the production of any heavy 
meson is described as the perturbative production of a colour-singlet QQ pair in a state 
with properly arranged the quantum numbers, according to the quarkonium state under 
consideration. For the production of J/%1) particles, the relevant partonic subprocess 
is 

7 + y ^ 35i[l]+^, (1) 

where we follow the standard spectroscopic notations, and the number in the brackets 
stands for the colour representation of the cc pair. The corresponding Feynman dia- 
grams are depicted in Fig. 1(a) (assuming also five possible permutations of the photon 
and gluon lines). The formation of a meson from the quark pair is a nonperturbative 
process. Within the nonrelativistic approximation which we are using, this probability 
reduces to a single parameter related to the meson wave function at the origin |7^(0)|^, 
which is known for the J / ip and T families from the measured leptonic decay widths. 

In addition to the above, we consider the colour-octet production scheme [8]-[ll], 
which implies that the heavy QQ pair is perturbatively created in a hard subprocess 
as an octet colour state and subsequently evolves into a physical quarkonium state 
via emitting soft (nonperturbative) gluons, which may be interpreted as a series of 
"classical" colour-dipole transitions. Although these transition probabilities can, in 
principle, be expressed in terms of field operators and therefore calculated, no such 
calculation exists to date. Thus, the transition probabilities remain free parameters, 
which are assumed to obey a definite hierarchy in powers of u, the relative velocity 
of the quarks in the bound system under study. This freedom is commonly used to 
estimate the colour-octet parameters by adjusting them to experimental data. 

In the case when the colour-octet QQ state is allowed, there appear additional 
contributions from the diagram of Fig. 1(a) and another set of diagrams shown in Figs. 
1(b) and 1(c) (including all possible permutations). The graphs shown in Figs. 1(a) 
and 1(b) correspond to the following partonic subprocesses 

7 + ^ ^ i5o[8]+^, (2) 
7 + ^ ^ '-5i[8] + ^, (3) 
7 + 5 ^ (4) 
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while the ones shown in Fig. 1(c) refer to the subprocesses 

7 + 5 ^ 'So[8i (5) 
7 + 5 ^ 'P,/[8]. (6) 

(The production of '^S'i[8] state in a 2 ^ 1 photon-gluon fusion process is forbidden by 
the colour and charge parity conservation.) Although the 2^2 subprocesses (2)- (4) 
are of formally subleading order in in comparison with (5)- (6), their role cannot 
be regarded as small. Since they contribute to very different regions of the phase 
space, they even can dominate over the 2—^1 subprocesses under the experimental 
conditions. We will discuss this in more detail in Sec. 4. Note also that the 2 — > 2 
subprocesses are indispensable for the inelastic events (i.e., the events with large final 
state hadron mass). 

The colour-octet matrix elements (usually denoted in the literature as (OlOslO)) 
responsible for the nonperturbative transitions in (2)-(6) are related to the fictituous 
colour-octet wave functions, that are used in calculations in place of the ordinary 
colour-singlet wave functions: 

(0|O8|0) = |^|7^8(0)r = |^47^|*8(0)|^ 

This equation applies to all ^-wave colour-octet states, and a similar relation holds 
for TZ'siO) and ^s{0) if the P-wave colour-octet states are involved. For the sake of 
uniformity, we will be consistently using the notation in terms of 7^(0) and TZ'{0) for 
both colour-singlet and colour-octet contributions. 

A generalization of the above formalism to the /cr-factorization approach implies 
two important steps. These are the introduction of noncoUinear gluon distributions 
(which we show here) and the modification of the gluon spin density matrix in the 
parton level matrix elements (which we explain in the next section). 

In the numerical analysis, we have tried two different sets of /cr-dependent gluon 
densities. In the approach of Ref. [41] based on a leading-order perturbative solution 
of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) [42] equation, the unintegrated gluon 
density J-'g{x, k"^, fi'^) is calculated as a convolution of the ordinary (collinear) gluon 
density G{x, /x^) with universal weight factors: 




Q{v,k'r,l^') = ^ Jo(2Va,ln(l/r?)ln(/.V4)), 4 < /.^ (8) 

Tj Krp 

where /j, is the factorization scale, Jq and Iq stand for Bessel functions (of real and 
imaginary arguments, respectively), and as = 3«s/7r. In the leading-order approxi- 
mation, the parameter ctg is connected to the Pomeron intercept a{0) = 1 + A, with 
A = 4 In 2. We use the value of A = 0.35 as it is accepted in many other our papers 
([43] and references therein). 

Another parametrization is based on the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) [44] evolution equation. This approach was originally proposed in [31], 
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[42] and is now frequently discussed in the literature [45], [46]. It recalls the kine- 
matic relation between the virtuality (f' and the transverse momentum /c^ of a parton: 
(f = k^/ {1 — x). Consequently, the ordinary gluon density G{x, q^) may be considered 
as giving the k"^ distribution also. In this approach, the unintegrated gluon density is 
derived from the "coUinear" density by simply differentiating it with respect to (f: 

Tg{x,k'^,fj,'^ = k'^) = ^G'(x,g2)|g2=fe2. (10) 

As the BFKL and DGLAP equations are known to collect different logarithms, we find 
it worth exploring the numerical consequences of this difference. For consistency, the 
same leading-order (LO) GRV set [47] was used in both cases as the input coUinear 
density. 

3 Matrix elements and differential cross section 

At first, we consider the relevant 2 — > 2 partonic subprocesses as they are given by 
the photon-gluon fusion mechanism. Let ki, k2, k^, pc, and Pc be the momenta of 
the incoming virtual photon and gluon, the outgoing gluon and the outgoing heavy 
(charmed) quark and antiquark, respectively; ei, €2, and 63 the photon and gluon 
polarization vectors; rric the quark mass, = 2mc, k=k2—k^, and a, 6, and c the 
eight-fold colour indices of the incoming gluon, the outgoing gluon, and the (coloured) 
cc state. We also introduce the projection operator J{S,L), which guarantees the 
proper spin and orbital angular momentum of the cc state under consideration. Then, 
the photon-gluon fusion matrix elements read 

xC^ tr{T''T^T^} [kl - 2{pcki)]-^[kl + 2(^5^3)]"^ + 5 permutations, (11) 
Mbiig ipg) = tr{ /i ( M + rrich^, J {S, L) } iG^^^ (/c2, £2, -/cs, £3, -k-, n) 

xC^ p'>dtr{T'T'^} [kl - 2{pcki)]-^[k^]-^ + 1 permutation. (12) 

In the above expression, G^^^ is related to the standard QCD three-gluon coupling 

G^'\p, A, q, k, v)^{q- pYg^^ + (k - qfg^^ + (p - kYg''^ (13) 

The factor represented by the SU{2i) generator matrix has to be replaced by the 
unity matrix if the outgoing cc state is a colour singlet. The coefficient stands for 
the normalization of the cc colour wave function and is equal to l/-\/3 and 1/2 for the 

singlet and octet states, respectively. 

The projection operator J{S, L) reads for the different spin and orbital angular 
momentum states [5]- [7]: 

J(%) = J(5=0,L=0) = 75(^ + ^0)/^;/', (14) 
J(35i) = J(5=1,L=0) = /(5,)(A + me)/m;/^ (15) 
j('Pj) = J{S=1, L=l) = i^s- me) 4{Sz) + me)/mf . (16) 

States with various projections of the spin momentum onto the z axis are represented 
by the polarization vector e{Sz). 
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In the nonrelativistic approximation which we are using, the relative momentum 
q of the quarks in the bound state is treated as a small quantity. So, it is useful to 
represent the quark momenta as follows: 

11 

Pc=2^V + ?' Pc=2^v--g, (17) 

where is the momentum of the final state quarkonium. The probability for the two 
quarks to form a meson depends on the bound state wave function ^'(g). Therefore, we 
multiply the matrix elements (11)-(12) by ^(?) and perform integration with respect 
to q. The integration is performed after expanding the integrand around q — 0: 

M{q) = M\g=o + {dM/dq'')\,=o + ... (18) 

Since the expressions for Ai\q=o, {dAi/dq°')\q=Q, etc., are no longer dependent on q, 
they may be factored outside the integral sign. A term-by-term integration of this 
series then yields [7] : 

*(g) = ^7^(x = 0), (19) 



(27r)3 ^ 

/(lf^«^) = -^^"(^^)^^'(^ = °)' (2°) 

etc., where IZix) is the spatial component of the wave function in the coordinate 
representation (the Fourier transform of '^[(q)). The first term contributes only to 5*- 
waves, but vanishes for P- waves because 7^p(0) = 0. On the contrary, the second term 
contributes only to P- waves, but vanishes for 5'-waves because T^s^) — 0- States with 
various projections of the orbital angular momentum onto the z axis are represented 
by the polarization vector e(L^. 

The polarization vectors t{Sz) and tiL^) are defined as explicit four- vectors. In 
the frame where the z axis is oriented along the quarkonium momentum vector, = 
(0, 0, Ip^I, £'^), these polarization vectors read: 

e(±l) = (1, ±1, 0, 0)7^2, e(0) = (0, 0, E^, KD/m^,. (21) 

States with definite 5*2 and L,, can be translated into states with definite total momen- 
tum J and its projection using the Clebsch-Gordan coefficients: 

e^\J,J,) = k^.L,-\,S,\J,J,)e^{L,Y{S,). (22) 

As far as the gluon spin density matrix is concerned, one should take into account 
that gluons generated in the parton evolution cascade do carry non-negligeable trans- 
verse momentum and are off mass shell. One can trace a clear analogy between the 
/cT-factorization approach and the equivalent photon approximation in QED showing 
that the polarization properties of the off-shell incoming gluon are similar to those of 
the incoming virtual photon. The off-shell photon spin density matrix is given by the 
full lepton tensor 

^ L"" = 8 p'^p'' - A{pk) (23) 

where p is the momentum of the beam particle, and k is the momentum of the emitted 
photon. A similar anzatz (Eq. (24), see below) is used in the /cr-factorization approach. 
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Neglecting the second term in the right hand side of (23) in the small x limit, k, 
one arrives at the spin structure e^^e*'' ~ p'^p^ . The latter can be rewritten in the form 

d^^kl^k^/lkrf, (24) 

where we have represented the 4-momentum k as k = xp + kx and applied a gauge 
shift ^ — k^/x. This formula converges to the usual X^e'^e*'^ = —g^" when 
kT — ^ 0. In the present calculations, we use Equ. (23) for virtual photons, and Equ. 
(24) for off-shell gluons. The expressions (23) and (24) merge with each other in the 
ultrahigh energy limit. The effect of the second term in (23) at HERA energies is 
found to be about 5 to 10 percent [43]. As we will see in Section 4, the presence of 
longitudinal components in the off-shell gluon spin density matrix has important impact 
on the quarkonium polarization. The final state gluon in (11)-(12) is assumed on-shell, 
J2€3^T ~ ~9^" ■ The evaluation of the traces in Eqs. (11)-(12) is straightforward and 
is done using the algebraic manipulation system FORM [48] . 

To calculate the cross section of a physical process we have to multiply the matrix 
elements squared by the gluon distribution functions and perform integration over 
the final state phase space. The multiparticle phase space Y{d?Pi/ {'2^Ei) S'^iJ^Pin ~ 
J2Pout) is parametrized in terms of transverse momenta, rapidities, and azimuthal 
angles: d^Pi/{2Ei)—(7r/2)dpf'pdyid(l)i/{27r). Let 0i and 02 be the azimuthal angles of 
the scattered electron and the initial gluon, and y^, y^, (f)^, and 03 the rapidities and 
the azimuthal angles of the J/ip particle and the coproduced gluon. Then, the fully 
differential cross section reads: 

«'"'e2|7^(0)|M ^ 1 ^ , ,,,2 

±U II X2i> ^ gpins O ^.^loura 

xTg{x2, k^T, /x^) dkfj, dklj. dpl^ dys dy^ ^ ^ (25) 

The phase space physical boundary is determined by the inequality [49] 

Gis,i,klklklml)<0, (26) 



where s = {ki + /C2) , t = {ki — p^p) , and G is the standard kinematic function 

The initial gluon momentum fractions Xi and X2 are calculated from the energy- 
momentum conservation laws in the light cone projections: 

(/i^i + /i:2)ii;+p|| = = m^Texp(|/^) + Ik^rlexp^y^), 

(27) 

{ki + /c2)b-p|| = X2y/s = m^T ex.p{-yi^) + |^3t| exp(-7/3), 

m^T = {iTi^ + Ip^tP)^^^- liere, we preserve exact kinematics and do not neglect the 
"small" light-cone component of the gluon momentum. The multidimensional integra- 
tion in Eq. (25) has been performed by means of the Monte-Carlo technique, using the 
routine VEGAS [50]. 

4 Numerical results and discussion 

We start the discussion by presenting a comparison between our theoretical calcu- 
lations and experimental data collected by the HI collaboration at HERA [30]. The 
collaboration reports on the measurement of a number of differential cross sections: 
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da/dp^,^, da/dp^^^, da/dz, da/dQ'^, da/dy, da/dy*, and da/dW, where pT,tp and p^^^ 
are the J/ip transverse momenta in the laboratory and 7*p center-of-mass systems, 
respectively, z is the J/ip inelasticity variable defined as 2; = {pxpPp) / {kiPp) , = —kf 
is the photon virtuality, y and y* are the J/ip rapidities in the laboratory and 7*p 
systems, and W is the j*p invariant energy. The data collected in the kinematic range 
2 GeV^ < Q2 < 100 GcV2, 50 GeV < W < 225 GeV, 0.3 < ^ < 0.9, p*T% > 1 GeV^ 
will be referred to as "sample 1" (see Fig. 2), while the data collected in the kinematic 
range 12 GeV^ < < 100 GeV^ 50 GeV< W <225 GeV, 0.3 < z < 0.9, p*T% > 1 
GeV^, p^.^ > 6.4 GeV^ will be referred to as "sample 2" (see Fig. 3). 

On the theoretical side, we have examined the sensitivity of model predictions to 
the choice of the gluon distribution functions, the renormalization scale in the strong 
coupling constant, the value of the quark mass, and the values of the nonperturbative 
colour-octet matrix elements. 

The effect of the different equations (BFKL versus DGLAP) which govern the 
evolution of gluon densities is found to be as large as a factor of 2 in the production 
cross section. This is illustrated by a comparison of dash-dotted and dashed histograms 
in Figs. 2 and 3, respectively. Note, however, that the parametrizations which we have 
used here represent two extreme cases, so that most of the other parametrizations 
available in the literature may be expected to lie in between our curves (see [35], [43]). 

A similar efi'ect is connected with the renormalization scale //.^ in the running 
coupling constant q;s(/x^). The calculations made with — k^j^ and /i^ — m^rp are 
represented by the dash-dotted and dotted histograms in Figs. 2 and 3. Note that 
the setting /^^ = k^j, is only possible in the /cT^-factorization approach, while it is 
meaningless in the coUinear calculations where the parton transverse momentum is 
neglected: /c2T = 0. 

The quark mass plays in the calculations two essentially different roles. The "cur- 
rent" mass rric is present in the expressions for the perturbative matrix elements (11)- 
(12). The "constituent" mass determines the phase space of the reaction via its connec- 
tion to the physical mass of the final state, = 2mc. However, it is worth pointing 
out that this connection is not strict. In fact, it may be violated by the effects of 
binding energy and internal quark motion [51]. The sensitivity of model predictions to 
the quark mass setting was examined in a toy calculation, where the "current" mass 
rric and the "constituent" mass m^/2 were treated as independent parameters. We 
found that the production cross section is only sensitive to the "constituent" mass but 
remains rather stable against variations in the "current" mass. In the rest of the paper 
we will be always using nic = m,^,/2 = 1.55 GeV. 

The contributions from the 2^1 colour-octet subprocesses are cut away by the 
experimental restriction z < 0.9. Turning to the 2 — > 2 colour-octet contributions, 
one has to take care about the infrared instability of the relevant matrix elements. 
In a rigorous approach, one has to consider the corresponding 2^1 subprocesses at 
next-to-leading order. Then, the interference between the LO and NLO contributions 
must cancel the divergent parts of the 2 — > 2 subprocesses. Such calculations have 
been performed in the collincar factorization in [52]. Since the corresponding results 
are not yet available in the fc^-facctorization, we use an approximate phenomenological 
approach. In order to restrict the 2 — > 2 subprocesses to the perturbative domain, we 
introduce the regularization parameter q^eg, so that all propagators are kept away from 
their poles by a distance not less than ?|eg- ^® argued that the nonperturbative 

parts of the 2 — > 2 subprocesses can be absorbed into 2^1 subprocesses; that is, when 
the emitted gluon is soft, one can consider the final state as represented by a single 
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particle rather than by two. This our suggestion is similar to the one made in papers 
[8]- [10], [53] in the usual factorization scheme. In this approach, the regularization 
parameter g^gg in the 2^2 processes and the nonperturbative colour-octet matrix 
elements in the 2^1 processes must be correlated [53] to avoid double counting 
between the hard and soft gluons in the final state (and so, to avoid sensitivity of 
the results to the choice of Q'feg)- The numerical results shown in Figs. 2 and 3 are 
obtained with setting Qxeg = 1 GeV^ (which may be regarded as the lower limit of the 
perturbative domain). 

The results shown in Figs. 2 and 3 are obtained with the nonperturbative colour- 
octet matrix elements of Rcf. [11]. If the values exractcd from the recent analysis [40] 
were used unstead, the contribution from the colour-octet states would be a factor of 
5 lower^. Irrespective to the particular choice of the nonperturbative matrix elements, 
the production of J/ ip particles at the conditions under study is not dominated by the 
colour-octet mechanism. This conclusion is also supported by an independent analysis 
[37] of recent ZEUS photoproduction data [54]. In general, the data are consistent 
with the predictions of the colour-singlet model and lie within the theoretical band 
provided by reasonable variations in model parameters. Although no need is seen in 
the colour-octet contributions, their presence does not lead to serious discrepancies in 
the visible kinematic 

We find it instructive to proceed the discussion by presenting a comparison with 
another data displayed by the HI collaboration in their earlier publication [55]. The 
data collected in the range 2 GeV^ < < gg GeV^ 40 GeV < < 180 GeV, z > 0.2 
constitute two correlated samples. In the "inelastic" sample, an additional cut on the 
final state hadron mass was applied, Mx > 10 GeV, while the sample without any 
additional cuts is referred to as "inclusive" . 

Although these data are, probably, less convincing in statistics, they are more useful 
and interesting from the theoretical point of view. Firstly, in the inclusive production 
sample, the high- 2; kinematic range 2 ~ 1 provides an access to the 2^1 colour-octet 
contributions. Secondly, the cut on the final state hadron mass Mx in the inelastic 
sample prevents the 2 — > 2 colour-octet matrix elements from infrared divergence. 

In Fig. 4, we present a comparison between the data and our predictions for the in- 
elastic production. The final state hadron mass Mx is calculated as the invariant mass 
of the proton remnant and the final state gluon produced in the 2 — > 2 subprocesses 
(l)-(4). The requirement that the final state mass Mx be large means that the final 
state gluon must be hard. Thus, the theoretical estimates of the colour octet contribu- 
tions are no longer dependent on the artificial regularization parameter g^eg- '^^^ 
see that the data are reasonably described within the colour-singlet production mech- 
anism, and the colour-octet contributions are not needed. Moreover, the colour-octet 
contributions to the high- 2; region look even superfluous. The overall situation seems 
to favour the low values of the nonperturbative matrix elements proposed in [40]. 

In the inclusive production sample (see Fig. 5), a significant deficit is seen in the 
^-distribution at 2; ~ 1, although the intermediate values of z are described quite 
well. In the J/ip transverse momentum spectrum, the deficit is only seen at low 

^The nonperturbative matrix elements given in Ref. [11] are as follows: |7?.35i(0)|^ = 8 x 10~^ 

GeV^, |7ei58(0)P = 8 X lO'^ GeV^, ITZsg^O)]'^ = 8 x lO'^ GeV^, \Tl'JO)\'^ = 7 x IQ-^ GeV^. 

1 

According to [40], the {^Sq) and (^Pj) matrix elements should be reduced by a factor of 5, while the 
{^Sf) matrix element should be reduced by at least a factor of 50 or enven set to 0. This violates 
the naive nonrelativistic QCD scaling rules, but is consistent with estimates obtained within the 
fcr-factorization approach by other athors [38] , [39] . 
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Pt.V") while there is no discrepancy at higher pT,ip values. These properties may be 
taken as an indication that the inclusive sample contains large contributions from 
diffractive processes. As a consequence, our calculations underestimate the absolute 
J/ip production rate by a factor of 3. 

We conclude the discussion by presenting our results on the J/ip spin alignment 
in Fig. 6. The data are expressed in terms of the Icptonic decay angular parameter a, 
which characterizes the azimuthal angle distribution measured in the J/ifj rest frame 
with respect to a given reference axis: dTu ~ 1 + q;cos(^). The cases a — 1 and 
a — —1 correspond to transverse and longitudinal polarizations of the J/ijj meson, 
respectively. Our calculations show that the fraction of longitudinally polarized mesons 
increases with increasing Q^, which is a consequence of the enhancement of the longi- 
tudinal component in the polarization vector of a virtual photon. Unfortunately, the 
experimental data are rather indcfinitive. In this case, we can derive no conclusions on 
the agreement or disagreement between the theory and experiment. 

In the situation which we consider, the effects of the gluon off-shellness are shad- 
owed by the effects of highly virtual photons and are not clearly visible. At a different 
situation, when either the initial photons are real or the process is mediated by gluons 
only, the J/ip polarization properties are found to be perfectly compatible with the 
predictions of the /cr-factorization approach, as it was demonstrated earlier in Refs. 
[36], [37] and [40] devoted to the analysis of J/ip photoproduction and hadroproduc- 
tion data. 

5 Conclusion 

Here we have addressed the issue of performing a global analysis of quarkonium 
electroproduction within the /c^-factorization approach. The state of the art has not 
yet reached the precise quantitative level. There are uncertainties connected with the 
choice of unintegrated gluon densities, the renormalization scale in the strong coupling 
constant, the inclusion of next-to-leading-order subprocesses, and the nonperturbative 
colour-octet transitions. 

At the same time, the /cr-factorization approach shows a number of important 
achievements. As a general feature, the model behavior is found to be perfectly com- 
patible with the available data on the production of various quarkonium states at 
modern colliders. The model succeeds in describing the p^-spectra of J/i/j, Xc and T 
mesons at the Fermilab Tevatron and provides a consistent picture of the production of 
J/ip mesons by real and virtual photons at HERA (as it was demonstrated earlier in 
Refs. [36], [37] and [40]). The model even succeeds in describing the polarization phe- 
nomena observed in both pp and ep interactions, thus providing an important insight 
for solving a long-term puzzle. 

As a result of the complex interplay of all theoretical uncertainties, the numerical 
analysis becomes rather umbiguous. The electroproduction data can be successfully 
described within the leading-order colour-singlet mechanism alone. No need is seen in 
the colour-octet contributions. At the same time, there are no serious contradictions 
with the data if the colour-octet contributions are included. The effects of the different 
unintegrated gluon distributions can hardly be separated from the effects connected 
with the running coupling constant. 
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Figure 1: Feynman graphs representing the photon-gluon fusion mechanism in the 
colour-singlet and colour-octet models. Only the perturbative skeleton is presented; 
the soft gluons corresponding to the nonperturbative colour-octet transitions are not 
shown. 

Figure 2: A comparison between the theoretical predictions and experimental data 
[30] in the kinematic range 1 < < 100 GeV^ 50 < < 225 GeV, 0.3 < ^ < 0.9, 
p*r^^ > 1 GeV^. Dash-dotted histogram, the colour-singlet contribution with BFKL 
gluon density and Q;s(A;|y); dashed histogram, the colour-singlet contribution with 
DGLAP gluon density and as{k2T)'i dotted histogram, the colour-singlet contribution 
with BFKL gluon density and Q;s(m|,^); solid histogram, the sum of the colour-singlet 
and colour-octet contributions, with BFKL gluon density, ^^(/clj.), colour-octet matrix 
elements as in [11], and ^^eg ~ 

1 GeV2. 

Figure 3: A comparison between the theoretical predictions and experimental data 
[30] in the kinematic range 12 < < 100 GeV^ 50 < W < 225 GeV, 0.3 < ^ < 
0-9) Pt% > 1 GeV^, p'^^^ > 6.4 GeV^. Dash-dotted histogram, the colour-singlet 
contribution with BFKL gluon density and ^^(/cly); dashed histogram, the colour- 
singlet contribution with DGLAP gluon density and q;s(/c|^); dotted histogram, the 
colour-singlet contribution with BFKL gluon density and Q;s(m|,^); solid histogram, 
the sum of the colour-singlet and colour-octet contributions, with BFKL gluon density, 
q;s(A;|j.), colour-octet matrix elements as in [11], and g^eg — 1 GeV^. 

Figure 4: A comparison between the theoretical predictions and experimental data 
[55] on the inelastic J/ip production in the range 2 < < 80 GeV^, 40 < < 180 
GeV, z > 0.2, Mx > 10 GeV. Dash-dotted histogram, the colour-singlet contribution 
with BFKL gluon density and as{k2T)', dashed histogram, the colour-singlet contri- 
bution with DGLAP gluon density and Q;s(A;|y); dotted histogram, the colour-octet 
contribution with BFKL gluon density, as{k2x) and colour-octet matrix elements as in 
[11]; solid histogram, the sum of the colour-singlet and colour-octet contributions. 

Figure 5: A comparison between the theoretical predictions and experimental data 
[55] on the inclusive J/ip production in the range 2 < < 80 GeV^, AO < W < 180 
GeV, z > 0.2. Dash-dotted histogram, the colour-singlet contribution with BFKL 
gluon density and «s(A;|j.); dashed histogram, the colour-singlet contribution with 
DGLAP gluon density and ^^(/clj^); dotted histogram, the colour-octet contribution 
with BFKL gluon density, as{klj.), colour-octet matrix elements as in [11], and q^Qg ~ 
1 GeV^; solid histogram, the sum of the colour-singlet and colour-octet contributions. 

Figure 6: A comparison between the theoretical predictions and experimental data on 
the J/ip spin alignment represented in terms of the decay lepton angular distributions. 
Left panel, kinematic range 1 (as in Fig. 2) [30], < 6.5 GeV^; right panel, kinematic 
range 1 (as in Fig. 2) [30], > 6.5 GeV^. 
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